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Dynamical dipole mode in heavy-ion collisions has been studied by using an isospin-dependent Boltzmann- 
Uehling-Uhlenbeck (IBUU) model. We investigate the dependence of centroid energy and strength of the y 
spectrum on beam energy, N/Z ratio and mass asymmetry. The calculated yield and angular distribution of the 
y-ray produced by dynamical dipole emission are consistent with the experimental data. The results show that 
the detailed study of dynamical dipole radiation can provide information on the isospin evolution of charge- 
asymmetric heavy-ion collisions around the normal nuclear density. 


Keywords: Heavy-ion collision, IBUU, Dynamical dipole mode 


DOI: 10.13538/j.1001-8042/nst.25.030501 


I. INTRODUCTION 


The giant dipole resonance (GDR), a collective mode of the 
nucleons, which plays an important role in nuclear physics re- 
search, has been widely studied experimentally [1—3] and the- 
oretically [4, 5] during the past decades. The pre-equilibrium 
GDR (dynamical dipole mode), generated by a collective mo- 
tion of protons against neutrons before the formation of a 
fully equilibrated compound nucleus (CN), has become a hot 
topic [6-17]. During the early stage of charge-asymmetric 
heavy-ion collisions, a large amplitude collective dipole oscil- 
lation is triggered along the symmetry axis of the strongly de- 
formed composite system; the oscillation emits prompt pho- 
tons, called dynamical dipole radiation, in addition to the pho- 
tons originating from a thermal excited GDR of the hot CN. 
The y yield from the pre-equilibrium GDR contains lots of 
information about the early stage of collisions when the CN 
is still in a highly deformed configuration. 

Several microscopic transport models, such as time- 
dependent Hartree-Fock (TDHF) [15], Boltzmann- 
Nordheim-Vlasov (BNV) [17], constrained molecular 
dynamics (CoMD) [16], and isospin-dependent quantum 
molecular dynamics (IQMD) [18—20], have been successfully 
applied to study properties of the dynamical dipole emission. 
A “bremsstrahlung” approach derived in Refs. [13, 14] was 
widely adopted to estimate the contribution of this pre- 
equilibrium collective dipole radiation emission. Here we 
use the isospin-dependent Boltzmann-Uehling-Uhlenbeck 
(IBUU) transport model [21] to make systematical calcu- 
lations of the y yield and the emission anisotropy of such 
prompt dipole radiation. 
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Il. DESCRIPTION OF THE METHOD 


The IBUU model [22] is widely used for treating heavy- 
ion collision dynamics and extracting important information 
of nuclear matter. The BUU equation reads 
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where do wy /d€) is in-medium nucleon-nucleon cross sec- 
tion, Vi» is relative velocity for the colliding nucleons. U is 
the mean field potential which includes the isospin-dependent 
term, f is the single particle distribution function, and p is the 
corresponding momentum. The isospin dependence is incor- 
porated into the model through the initialization and nuclear 
mean field. The mean field including isospin symmetry term 
is parameterized as 


U(p,72) = A(Z) HAE Vif? (p5), — Q) 
Po po 

where po is the saturation nuclear matter density, and p = 
Pn + Pp is density of nucleon; pn and p, are the neutron and 
proton densities, respectively. The coefficients a, b and c are 
parameters for the equation of state (EOS). Here we choose 
the parameters corresponding to the soft EOS with the com- 
pressibility K of 200 MeV (a = —356 MeV, b = 303 MeV, 
and o = 7/6). The form of single particle symmetry poten- 
tial [23-25] reads 
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where Cs is the symmetry energy coefficient depending 
on the different densities of neutrons and protons in nuclear 
medium. In this paper, we use the measured symmetry en- 
ergy of C's — 35.19 MeV and, also, a soft symmetry energy 
(y — 0.5), which was found to describe the experiment data 
reasonably well in Ref. [17, 26]. As usual, the test particle 
method is applied in the IBUU simulations and the in-medium 
nucleon-nucleon cross section is taken into account [27]. For 
calculating dynamic dipole emission, the emitted nucleons 
with their local density below po /8 are excluded. 

During pre-equilibrium dipole oscillation, a larger initial 
dipole moment triggers higher amplitude isovector oscilla- 
tions, which increases the chance of a clear experimental 
observation. The giant dipole moment in coordinator space 
D(t) is defined as [14] 
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where A= N + Z, N = Np + Ni, Z = Zp + Z are the 
total number of participating nucleons, neutrons, and protons 
from the projectile (p) and the target (t), respectively; X (t) is 
the distance between the centers of protons and neutrons; and 
Rp and R, are the total coordinates of center-of-mass system 
for the protons and neutrons, respectively. And the canonical 
conjugate momentum of D(t) is written as 
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where P, and P are the total momenta of center-of-mass 
system for the protons and neutrons, respectively. The initial 
(t = 0: touching configuration) giant dipole moment can be 
expressed as 
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where Rp and R; are the radii of the projectile and target, 
respectively [13]. 

Then, as given by the bremsstrahlung approach, the » yield 
can be extracted as [13, 14] 
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where E, = hw is photon energy and | D”(w) |? is the 
Fourier transform of the dipole “acceleration” D"(uw) = 
fS m D" (t)e***dt. For each event, to represents the onset- 
time of the collective dipole response (phase-space spiraling) 
and tmax is the “damping time". 

Recently, a clear anisotropy for the angular distribution of 
the total ^; spectrum has been observed [10, 11]. For the 
dipole oscillation just along the beam axis, the angular dis- 
tribution is expected to be given by the Legendre polyno- 
mial expansion M (0) ~ sin? 0 ~ 1 + aa P»(cos 0) (a2 is the 
anisotropy parameter and 0 is the polar angle between the 
photon-emitting direction and oscillation axis). For the dy- 
namical dipole mode, the prompt dipole axis rotates during 
the emission. 
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First, in a very small time interval At, p; and vy are de- 
fined as the initial and final angles of the oscillation axis with 
respect to the beam axis, respectively. So Ay = yy — yi is 
the rotation angle during the time interval At. From preced- 
ing text, the angular distribution in At can be averaged over 
the angle Ay as M (0) ~ 1 — (1/4 + 3x/4) P2(cos 0), where 
x = cos(o poi) sin(pr — vi)/(or — vi) [17]. Then, for the 
whole time of the dynamical dipole, we can extract the angu- 
lar distribution by using a weighted form for every time step 
and summing over all of them: M(0) = Y, 8;M (0, ®;), 
where 6; is the rotation angle and 5; is the radiation emis- 
sion probability. From Eq. (7) we have 8; = P(t;) — P(t; 1) 
where P(t;) — p D"(t) |? dt/ B, and Pot is given by 
P(tmax), which is the total emission probability at the final 
dynamical dipole damped time. Finally, the anisotropy pa- 
rameter is: 
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The angular distribution of the prompt y emission has been 
extracted by using the formalism from Ref.[17]. 


IIl. RESULTS AND DISCUSSION 
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Fig. 1. (Color online) The results in 16 MeV/nucleon 36 Ar + 967; 
with b = 4 fm: time evolution of dipole moment D(t) in coordinate 
space (a), dipole phase-space correlation (b), and comparison of the 
"y yield (c) and angular distribution (d) between our calculations and 
experimental data under the same conditions. 


Figures 1(a) and 1(b) show the time evolution of the dipole 
moment D(t) and the phase space correlation (spiraling) 
between D(t) and DK(t) for the 16 MeV/nucleon ?9Ar + 
967r with b = 4fm. One sees that the large amplitude of 
the first oscillation decaying rapidly within several periods, 
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Fig. 2. (Color online) Evolution of (a) the total strength (P), (b) the 
width (W) and (c) the central energy (Ee) of y spectrum in the 36 Ar 
+ 87r reactions at Eveam = 12, 16, 20 and 24 MeV/nucleon with 
b = 0, 2 and 4fm. 


and the whole process continues for about 300 fm/c. Espe- 
cially, Fig. 1(b) shows the collective oscillation initiated at 
the very early time, which corresponds to the touching con- 
figuration value of the collision. So, in this paper, we de- 
fine tmax — to = 300fm/c to make sure that the collective 
oscillation is basically over. In Fig. 1(c), the calculated pho- 
ton yields for the 16 MeV/nucleon ?9Ar + 997r are compared 
with the experimental data [10]. Fig. 1(d) shows the angular 
distribution, both the simulation and experimental results are 
averaged from data obtained at impact parameters of b — 0, 
2 and 4 fm. Our calculation results can reproduce the experi- 
mental data fairly well, indicating that the model is successful 
in describing the y radiation of dynamical dipole mode. 


So far, both experimental data and theoretical calculations 
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Fig. 3. (Color online) Anisotropy parameter of ?9 Ar + °°Zr at dif- 
ferent incident energies and impact parameters. 


have shown that y spectrum of the dynamical dipole mode 
is dependent on beam energy of the reaction system. In this 
paper, detailed results are drawn from the dynamical emission 
of 12, 16, 20 and 24 MeV/nucleon 3° Ar + 967r with different 
impact parameters. Using the Lorentz fitting, we extract the 
central energy (E<), the width (W), and the total strength (P) 
of the y spectrum (Fig. 2). 

As shown in Fig. 2(a), the strength has a maximum value 
around Eteam = 16 MeV/nucleon, which is the same as the 
results in Ref. [14]. At low incident energies, the dipole os- 
cillation increases fiercely with E545; while at high incident 
energies, the damping properties of the dipole moment be- 
come stronger and fast nucleon emission takes place which 
make the dipole oscillation become weaker [26]. The damp- 
ing properties and the nucleon emission also lead to the evolu- 
tionary trend of dynamical dipole emission width (Fig. 2(b)). 
In Fig. 2(c), the dynamical vibration has a high frequency 
at low incident energies. This means a faster isospin equili- 
bration in the collisions at lower beam energies, which was 
concluded similarly in Ref. [28]. Also, the strength of y radi- 
ation increases obviously with the impact parameter, because 
of the increasing dipole moment on the direction perpendic- 
ular to the beam. The central energy of y spectrum has the 
same behavior, because part of the beam energy transfers to 
the rotational energy of CN when the reaction happens away 
from the central collision. 

Figure 3 shows the calculation results of anisotropy param- 
eters (a2) of ?6Ar + °°Zr at different incident energies and 
impact parameters. The a» value decreases faster at higher 
incident energies, because higher beam energy causes faster 
rotation of the oscillation axis. The same thing happens as the 
impact parameter becomes larger. 

The dynamical dipole emission originates from different 
N/Z ratios (Ry /z) between target and projectile, which leads 
to the process of charge equilibration. Fig. 4 shows the E and 
P of the y spectrum produced by reaction systems of ?9Ar 
O67, (Ry jz = 1.400), 36 Ar + °°Mo (Ry jz = 1.286) and 
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Fig. 4. (Color online) Calculated (a) total strength (P) and (b) cen- 
tral energy (Ec) of the y spectrum produced by the reactions of 
36 Ar + 997r (Ryyz = 1.400),  Ar+°®Mo (Ry /z = 1.286) and 
36 Ar+°°Ru (Ry/z = 1.182) at 16 MeV/nucleon with different im- 
pact parameters. 
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Fig. 5. (Color online) Calculated (a) total strength (P) and (b) cen- 
tral energy (E) of the y spectrum produced by the reactions of ?9 Ar 
+ ©°Rb (6mass = 2.472), 8 Ar + °°Zr (ômass = 2.667) and ĉĉ Ar + 
losg (mass = 2.861) at 16 MeV/nucleon with different impact 
parameters. 
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6 Ar + ?9Ru (Ry jz = 1.182) at 16 MeV/nucleon. Fig. 4(a) 
demonstrates a sharp monotonic increasing behavior in the 
correlation between Ry jz and P. The large difference of 
dipole moments at the incident channel is responsible for this 
behavior. In Fig. 4(b), the Ee varies just a little at different 
values of Ry z, but increases with the impact parameters. 
Besides, at each impact parameter, the anisotropy parameters 
(a2) of different reactions are almost the same (not shown). 

To study the influence of mass asymmetry, the reac- 
tions of ?9Ar + Rb, 96Ar + 967r and 36Arp + Tce at 
16 MeV/nucleon are calculated. The N/Z ratios of these 
systems are the same (Ey,z = 1.400), but the parameters 
of mass asymmetry (mass = MASSaarget "i MASS projectile) are 
2.472, 2.667 and 2.861, respectively. Compared with the re- 
sults in Fig. 4, an opposite behavior is seen in Fig. 5: the vari- 
ation of P is small and the E. has a decreasing tendency with 
a larger dass. But all of them increase monotonously with 
the impact parameter. The yields of dynamical dipole radi- 
ation are mainly dependent on incident energy and the N/Z 
ratio. And the results of anisotropy parameters at different 
values of dqass are similar to the case of charge asymmetry. 
In other words, the major effects on the angular distribution 
are the beam energy and the impact parameter. 


IV. SUMMARY 


In this paper, we use the IBUU model to study the 
properties of the dynamical dipole mode via a collective 
bremsstrahlung mechanism in fusion reactions with charge- 
asymmetric systems. We compare the y spectrum and angu- 
lar distribution of the prompt photon radiation with the exper- 
imental data, and show the features of this dynamical dipole 
emission with low energy reactions of different systems. We 
study the total emission probability, centroid energy and an- 
gular distribution of the pre-equilibrium dipole oscillation and 
their dependence on incident energy, N/Z ratio and mass 
asymmetry. Our results demonstrate that the total emission 
probability of the dynamical dipole mode has a strong depen- 
dence on beam energy and charge asymmetry in the entrance 
channel, besides on mass asymmetry. The central energy of 
this pre-equilibrium dipole oscillation is more sensitive to the 
beam energy and mass asymmetry than the case of charge 
asymmetry. And just the beam energy has an influence on the 
anisotropy of the prompt photon radiation. But, almost, all of 
the features are sensitive to the impact parameters. So, this 
work indicates that detailed study of dynamical dipole radia- 
tion can be used as a probe to extract information on the early 
stage of charge-asymmetric heavy-ion collisions around the 
normal nuclear density. 
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